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Abstract 

Using the simple relation between the np scattering and bound state 
wave functions, the spin-triplet contributions to the pp — ► 7r + (pn) and 
backward dp — ► p(pn) reactions are estimated for low np excitation en- 
ergies. The good agreement with the pion production data at 400 and 
450 MeV over a range of angles shows that spin-singlet final states can- 
not be more than of the order of 10%. On the other hand it is seen that 
the spin-singlet states must provide about 30% of the strength in pion 
production at 1 GeV as well as deuteron-proton backward scattering at 
1.6 GeV. 
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It is well known that when a scattering wave function is extrapolated to a 
bound state pole then the result is proportional to the bound state wave function. 
Less appreciated is the fact that the constant of proportionality is independent of 
the form of the potential |l], 0. Thus for single- channel S-wave neutron-proton 
scattering, 

Jto {-^*M)e«tf-»(r)}= W r). (1) 

Here at = y/rne = 0.232 fin" 1 , where m is the average nucleon mass, e the 
deuteron binding energy, and 5 the uncoupled S-wave triplet phase shift at relative 
momentum k. Though this relation is only true at the deuteron pole, nevertheless 
for cm. energies below about 40 MeV and internucleon distances r less than 
about 1.5 fm the 3 Si solution to the Paris neutron-proton potential p| empirically 
behaves like 

to an accuracy of typically 5% ||. Thus for large momentum transfer reactions 
such as 7r [||] or i] |1] production, which are sensitive to the values of the wave 
functions at short distances, one can get a simple estimate of the amplitude for 
the production of a particle c, with a neutron-proton pair in the final state in 
terms of the amplitude for deuteron (d) production, viz. 



/ 27T7TZ 

M(ab -> c{np} k ) « -J — e~ iS M{ab -> cd) . (3) 

V a t {k z + af) 

Though there is no bound state amplitude to set the scale, the momentum de- 



pendence of the S-wave Paris singlet wave function is given by a similar ^jk 2 + a 2 s 
factor to that of eq.(l), but with a different pole parameter a s = — 0.104 fm -1 
[El p[ . This factor in turn determines the energy variation of the production of np 



singlet final states. Such a simple model for triplet and singlet production gives 
a good description of the total cross sections for np — > nprj |IJ and pp — > pmi + 
and pd — > pdn @ near threshold. It is the aim of the present work to show that 
the model is equally successful at higher energies providing the neutron-proton 
system is constrained to be at low excitation energy. 

Choosing the same normalisation for bosons and fermions, the cm. differen- 
tial cross section for the production of an N-body final state is 

Here pj and Ej are the momentum and energy of the final-state particle j and 
Pi the incident momentum corresponding to a total cm. energy -^i, whereas pi 
and pf are the total initial and final four- momentum vectors. 

By introducing the ansatz of eq.(|3|) into the phase space and making kinematic 
approximations which are very good for high beam energies and low excitation 
energies Q = k 2 /m in the final np system, the spin-triplet contribution to the 
three-body cross section is predicted to be 

rfV <ab - c{np}) = 4^ f , ^(ab - cd) . (5) 



dtt c dx y L J/ P(-1) 2tt(1 + x) dtt c 

The dependence on the energy of particle c is mainly through the explicit 
function of x, where 

x=^- = - 2 - (6) 

p{x) is the cm. momentum of particle c in the c{np} system, becoming equal 
to that in the cd system when x = —1. The ratio p(x)/p(—l) is independent of 
angle and depends weakly upon x and the incident beam energy. 
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The most extensive measurement of the pp — > ir + pn cross section and analysing 
power was carried out several years ago at TRIUMF at 400 and 450 MeV at sev- 
eral fixed laboratory angles Unfortunately the data are no longer available in 
numerical form but it is possible to read values from the published work f| and 
the original hand-drawn graphs || to a fair accuracy. It is known over a wide 
energy domain from threshold up to the A-region that the final spin-triplet states 
dominate the reaction || and indeed it was noted that the pp — > 7i + pn analysing 
powers follow closely those of pp — > ir + d at the same cm. pion production angle 

On the basis of eq.©, the triplet contribution to the cross section ratio is 

A crucial point is that the same pp — > ir + d cross sections employed in the nor- 
malisation in the original data analysis has to be used when constructing the 
experimental value of the ratio. This information is indeed available ||. 

In fig. 1 all the data at the two energies and eleven laboratory angles are 
plotted against Q and their mutual consistency and the agreement with the pre- 
dictions of eq.(0) at small Q is impressive. The immediate conclusions that one 
can infer from this comparison are: 

1. Since any contributions from spin-singlet and P-wave final states to the 
semi-inclusive cross sections are incoherent and non-negative, the exper- 
imental points should not lie below the theoretical curves. Violations of 
this bound occur at at large laboratory angles where, due to kinematic 
effects, the experimental resolution is poorer. The isolation of unbound 
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final neutron-proton states from the deuteron peak, which is not shown but 
which is used to normalise the cross section, becomes more difficult under 
such conditions. 

2. Taking into account that eq.(2) may be in error by of the order of 5%, the 
possible amount of spin-singlet final np states must be generally less than 
about 10% of the spin-triplet. 

3. There is evidence from the more reliable small angle points that, above an 
excitation energy of 10 MeV, the curves represent a smaller fraction of the 
data, as would be expected from the neglect of P and higher waves in the 
np system. 

4. It has been pointed out Q that the angular distribution of the proton 
analysing power follows closely that of pp — > ix + d at low excitation energy. 
Apart from the larger angle points, the agreement with eq.(^) is good over 
a range of laboratory angles, and so this provides evidence for a similar 
behaviour for the differential cross section. However the resolution problem 
must be clarified to make this more quantitative. Fortunately the approach 
will soon be tested in far greater detail through the new high accuracy 
pp — > ir + pn TRIUMF data currently being analysed [f|. 

Turning now to a second application of the technique, the momentum transfer 
in elastic deuteron-proton scattering in the backward direction increases very fast 
with beam energy. For deuteron laboratory kinetic energies below Td « 200 MeV 
the process is dominated by one-neutron-exchange and the cross section then 
depends on the deuteron wave function at high Fermi momenta. On the other 
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hand, for Td > 800 MeV the process is expected to be driven by pionic degrees of 
freedom and virtual A-excitation JFJ. The gap region between these two models 
is less well understood. 

The dp — > pd cross section and deuteron tensor analysing power T 20 have 
been measured at 180° for 300 < Td < 2300 MeV and data are also available on 
the corresponding three-body reaction dp — > p(pn) at 400 and 1600 MeV f|. At 
1600 MeV the momentum transfer is so large that, just as for pion production, the 
reaction should be dominated by short-range physics. Since the energy resolution 
in the data is not as good as for the pion production experiment previously 
discussed, smearing of the predictions for the final deuteron and proton-neutron 
pairs over the experimental resolution function is required. Using eq.([5]), the 
summed cm. cross section can be written as 

da {dp -»• p{pn}) = + x) + ^ T^ ( x )| ^{dp ^ pd) dVt p dx , 

(8) 

where the delta-function reflects the two-body pd final state and the term involv- 
ing the theta- function the three-body triplet p(pn). 

Convoluting this distribution with an experimental resolution function in Q, 
taken to be Gaussian with o = 2.35 MeV, leads to the spin-triplet prediction 
shown in fig. 2a for Td = 1600 MeV. It must be stressed that it is the area of 
the elastic deuteron peak which fixes the normalisation of the triplet break-up 
contribution through eq.(|8|). 

Unlike the low energy pion production data shown in fig. 1, our estimate 
leaves some space for final np singlet states. The shape of such a cross section is 
dominated by a similar ^fxj[x + a^/af) form but, because of the smaller value 
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of a s , it is more peaked in the near-threshold region. Tolerable agreement with 
the data is achieved with a singlet final-state cross section about 0.4 ± 0.05 of 
the triplet at large Q, which is essentially consistent with a purely statistical 
population. 

In contrast, at = 400 MeV it is not possible to get a good fit to the data 
shown in fig. 2b even after including singlet np final states. The enhancement at 
small Q is incompatible with the data at larger values since np P-waves can only 
increase the cross section there. This breakdown of the model is not unexpected 
since at 400 MeV the one-neutron-exchange term is still the largest contribution 
and the short-range limit employed here is of limited value. 

The difference in behaviour at the two deuteron beam energies is also seen in 
the deuteron tensor analysing power T 20 . Values of this quantity were measured 
for the sum of the first six bins in np excitation energy as well as for elastic proton- 
deuteron scattering f8fl. In the 1600 MeV region the difference AT 20 between these 
two seems to be energy independent and consistent with zero. On the other hand 
both the elastic T 2 o and AT 2 o vary fast with energy in the 300 - 600 MeV region, 
indicating a dependence on the details of the deuteron wave function. It is clearly 
not the high momentum transfer domain where our approach can be used in a 
naive way. 

The most plausible model for large angle elastic pd scattering at high ener- 
gies involves the virtual production and absorption of pions J7[ and the striking 
difference between the singlet /triplet ratio determined from pion production at 
400/450 MeV and the dp — > p(pn) reaction at an equivalent proton beam energy 
of 800 MeV could bring this into question. There is however evidence from pion 
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absorption on 3 He that above the A resonance the singlet /triplet ratio is much 
larger than at lower energies 0. This is backed up by the pp — ► 7i + pn data 
taken at 1 GeV || and shown in fig. 3. The absolute resolution is here a little 
poorer than for the 1.6 GeV dp — » p(pn) experiment |J but, when the latter data 
are degraded slightly by including small admixtures from neighbouring bins, the 
shapes of the two data sets are indistinguishable. Kinematic differences due to 
the variation of the pion/proton momentum with Q are negligible at such high 
energies. Also shown are the fitted curves of fig. 2a with the slightly increased 
smearing width (a = 2.60 MeV). Pion production at 1 GeV is therefore also 
consistent with a singlet/triplet ratio of around 0.4 ± 0.05. 

This steady increase of the singlet /triplet ratio with beam energy at the A- 
resonance and above is confirmed in the analysis of p(p, tt + )X results at 592 MeV 
and 9 = 0° [K|. A good fit to this data is achieved with a ratio of 0.25 ± 0.05. 

We have shown that a plausible description can be given for high momentum 
transfer reactions leading to unbound np scattering states at low excitation en- 
ergy Q in terms of the corresponding cross section for producing a deuteron. In 
the limit that the transition operator is of short range our evaluation is weakly 
model-dependent. Our principal conclusions are that any spin-singlet final states 
in pp — > 7T + pn at 400 and 450 MeV must be 10% or less of the triplet, though 
this steadily increases to 20% and 40% at 600 and 1000 MeV respectively. This 
final figure, which is essentially consistent with a statistical spin factor of 3:1, 
also reproduces well the 1600 MeV dp — > ppn data, giving some support to the 
idea that this reaction is driven by virtual pion production. To go further would 
require the introduction of an explicit reaction model as well as a fuller evaluation 
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of the final neutron-proton wave function including the deuteron D-state. Our 
calculation has demonstrated that such a refined model should give reasonable 
answers. A systematic experimental study of the pp —>■ tc + X reaction in the for- 
ward direction as a function of energy should show clearly the evolution of the 
singlet /triplet ratio and the importance of non-A contributions to the reaction. 
The value of such experiments would be enhanced if they could be carried out 
with polarised beam and target since, by angular momentum conservation, in 
the forward direction an initial longitudinal spin-spin correlation is propagated 
to the final neutron-proton system. 

Considerable assistance in using the data of ref.|| was given by W.R. Falk and 
this made our work much easier. The values of the nucleon-nucleon S-wave singlet 
and triplet scattering state wave functions were kindly furnished by B. Loiseau 
and we have benefited greatly from his advice in their interpretation. Valuable 
discussions with V. Koptev are also gratefully acknowledged. This work has been 
made possible by the continued financial support of the Swedish Royal Academy 
of Science and one of the authors (CW) would like to thank them, the The 
Svedberg Laboratory and the Kernforschungzentrum Jiilich for their generous 
hospitality. 
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Figure 1: The ratio R of the cm. differential cross section for pp — > n + pn 
divided by that for pp — > n + d at the same pion cm. angle, as defined by eq.(^). 
The TRIUMF experimental data [|J] were taken at 400 MeV (open symbols) and 
450 MeV (closed symbols) at a total of 11 laboratory angles. The predictions of 
eq.(0) at these two energies are shown as the broken and solid line respectively, 
the only difference between them being due to the change in the variation of the 
momentum factor. 
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Figure 2a: Laboratory differential cross section for the p(d,p)pn reaction in the 
backward cm. direction at Td = 1600 MeV as a function of the excitation 
energy Q in the unobserved final np system ||. The binning is such that 1 
channel corresponds to 2.69 MeV. The form of eq.([8D was smeared over a Gaussian 
resolution function with a = 2.35 MeV and a binning correction applied before 
being compared to the data. The broken line is the fitted deuteron elastic peak 
and the dot-dashed the resulting prediction for the triplet final np state. The 
amount of spin-singlet contribution required to give the overall fit (solid line) is 
shown dotted line. 
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Figure 2b: Laboratory differential cross section for the p(d,p)pn reaction in the 
backward cm. direction at Td = 400 MeV as a function of the excitation energy 
Q in the unobserved final np system ||. The binning is such that 1 channel 
corresponds to 0.89 MeV. The form of eq.(|8|) was smeared over a Gaussian reso- 
lution function with a = 0.52 MeV and a binning correction applied before being 
compared to the data. The broken line is the fitted deuteron elastic peak and the 
dot-dashed the resulting prediction for the triplet final np state. No satisfactory 
overall fit to both the low and high excitation energy data (solid line) can here 
be achieved by adjusting the proportion of spin-singlet final states (dotted line). 
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Figure 3: Laboratory differential cross section in arbitrary units for the p(p, ix + )X 
reaction at 9 = 0° and T p — 1.0 GeV (stars) as a function of the excitation en- 
ergy Q in the unobserved final np system ||. These are compared with the 
1.6 GeV p(d,p)X data of ref. || (circles) whose energy resolution has been de- 
graded slightly by replacing the cross section a(n) in the n'th bin by 0.15 a(n + 
1) + 0.7a(n) + 0.15 a{n — 1). The fitted curves are identical to those in fig. 2a 
but with a smearing parameter of a = 2.60 MeV. 
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